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ABSTRACT 

The  translational  and  rotational  motions  of  two  prolate  spheroids 
sedimenting  in  a  viscous  fluid  have  been  determined  by  the  method  of 
reflections.  No  restrictions  are  imposed  on  the  spheroid  orientations  or 
relative  sizes.  As  is  the  case  in  many  mobility  problems,  the  method 
converged  rapidly  for  all  but  almost  touching  configurations.  The  results 
extend  earlier  work  on  special  cases  such  as  Wakiya's  work  on  horizontal 
orientations  and  agree  with  Gluckman  et.  al.  and  Liao  and  Krueger's  boundary 
collocation  solution  of  axisymmetric  problems.  Analysis  of  sedimentation  with 
inclined  axes  and  mirror  symmetric  geometry  reveal  both  periodic  and  single- 
encounter  particle  trajectories.  The  calculation  of  the  Beparatrix  between 
the  two  behaviors  required  the  use  of  the  higher  reflections  introduced  in 
this  work. 

AMS  (MOS)  Subject  Classifications:  76D05,  35Q10 

Key  Words:  Sedimentation,  Spheroids,  low-Reynolds-number ,  Hydrodynamic 
Interaction. 

work  Unit  Number  2  -  Physical  Mathematics 


SIGNIFICANCE  AND  EXPLANATION 


The  calculation  of  hydrodynamic  interactions  between  particles  is  needed 
for  the  understanding  and  control  of  many  natural  and  manufacturing  processes, 
for  instance,  those  involving  sedimentation,  colloidal  stability  or  suspension 
rheology.  In  these  applications,  the  external  forces,  torques  and  the  ambient 
velocities  are  known  a  priori  and  the  problem  is  to  calculate  the  resulting 
translational  and  rotational  motions.  In  practice,  since  the  governing 
differential  equations  require  knowledge  of  these  motions  for  the  boundary 
condition,  one  has  to  solve  first  for  the  forces,  torques  and  dipole  moments 
in  a  collection  of  translational  and  rotational  problems  and  then  invert  them 
to  obtain  the  desired  motions. 

The  purpose  of  this  report  is  to  show  that  these  problems  can  be  solved 
directly.  Explicit  calculations  and  comparisons  with  other  techniques  are 
shown  for  the  sedimentation  of  two  spheroids  in  a  viscous  fluid.  Deviations 
from  the  settling  behavior  of  spherical  particles,  such  as  particle  drift,  are 
highlighted. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


SEDIMENTATION  07  TWO  ARBITRARILY  ORIENTED 
SPHEROIDS  IN  A  VISCOUS  FLUID 

Sangtae  Kia 

1.  INTRODUCTION 

Suspensions  of  prolata  spheroids  have  played  an  important  role  In  the 
theoretical  development  of  suspension  rheology.  Suoh  suspensions 
exhibit  non-Newtonian  behavior  through  the  interaction  between  the  flow 
field  and  Brownian  motion  (Glesekus  1962,  Brenner  1972,  Hinoh  and  Leal 
1972.  However,  rigorous  derivation  of  the  material  funotions  to  date 
have  been  restrloted  to  the  dilute  limit,  partly  because  of  the  lade  of 
information  on  multl-partiole  hydrodynamio  interactions.  Existing 
information  on  partlole-particle  interactions  is  limited  to  certain 
geometries  at  large  parti ole-partlole  separations  (Vakiya  1965)  or 
Special  configurations  (Gluckman  at.  al.  1971;  Liao  and  Krueger  1980). 

(This  report  continues  the  work  in  (Kim  1984),  MRC  Technical  Summary  Report 
*2643.) 

New  results  are  presented  here  which  describe  the  interactions 
between  two  spheroids  with  arbitrary  configurations  and  all  but 
almost-touohing  separations.  Explicit  examples  are  worked  out  to 
illustrate  phenomena,  such  as  the  evolution  of  particle  geometry,  which 
are  not  found  in  the  corresponding  problem  for  spheres.  The 
computational  technique  which  is  based  on  the  method  of  reflections 
(Happel  and  Brenner  1965,  Felderhof  1977),  was  found  to  eonverge  rapidly 
for  the  sedimentation  and  related  mobility  problems.  The  improved 
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results  for  the  mobility  functions  were  essential  in  accurate 
calculation  of  particle  trajectories,  especially  in  the  near-field 
Interactions. 

In  Section  2,  the  techniques  for  calculating  hydrodynamic 

interactions  which  were  developed  in  an  earlier  note  (Kim  198H)  are  used 

to  recover  Wakiya’s  (1965)  results  for  the  resistance  problem.  In 

Section  3,  sedimentation  and  angular  velocities  are  calculated  to  0(R f^) 
“8 

and  0(R  )  respectively,  where  R  is  the  center  to  center  separation 

between  the  spheroids.  An  advantage  of  the  present  method  is  that  it 
bypasses  the  usual  procedure  of  calculation  and  inversion  of  the 
resistance  problem.  As  outlined  in  Kim  (1983),  the  sedimentation 
problem  is  solved  directly,  without  solving  a  collection  of  subsidiary 
problems  on  translating  and  rotating  spheroids.  Problems  solved  include 
sedimentation  along  and  perpendicular  to  the  line  of  centers  and  the 
evolution  of  configurations  for  spheroids  with  Inclined  axes. 

2.  HYDRODYNAMIC  INTERACTION  BETWEEN  TWO  STATIONARY  SPHEROIDS 
In  this  section,  the  method  will  be  used  to  recover  Wakiya’s  (1965) 
calculations  for  the  drag  on  a  spheroid.  Figure  1  shows  the  geometry 
used  by  Wakiya  (1965).  In  order  to  simplify  the  final  expressions,  he 
restricted  his  analysis  to  two  identical  spheroids  with  both  axes  placed 
horizontally  (with  gravity  acting  in  the  negative  z  direction).  The 
drag  and  torque  on  spheroid  1  was  calculated  for  the  case  where  both 
spheroids  were  translating  (without  rotating)  in  the  negative  x 
direction  in  a  quiescent  fluid. 

In  the  terminology  of  the  general  literature,  this  is  called  a 
resistance  problem.  The  translational  and  rotational  velocities  are 
specified  and  the  drag  and  torques  are  to  be  found.  The  mobility 


problems  pose  the  inverse  question,  i.e.,  forces  and  torques  on  the 
particle  are  specified  and  the  translational  and  rotational  motions  are 
to  be  determined.  The  latter  problem  occurs  more  frequently  in  the 
modeling  efforts  of  diverse  fields.  Sedimentation  and  diffusion 
problems  in  suspension  rheology  and  hydrodynamic  interactions  in  the 
Rouse-Zimm  theory  all  require  the  solution  of  a  mobility  problem. 
Specific  applications  can  be  found  in  the  following  samples  from  an 
extensive  list:  Glendinning  and  Russel  (1982),  Batchelor  (1976),  and 
Bird,  Hassager,  Armstrong  and  Curtiss  (1977). 


The  sedimentation  problem  (e  nobility  problem)  and  the  problem  of 
calculating  the  drag  on  atatlonary  objects  (a  resistance  problem)  are 
reciprocals  of  each  other  up  to  0(R-^) .  This  simple  situation  does  not 
hold  at  higher  orders  because  torques  are  present  in  the  reslstanoe  but 
absent  in  the  sedimentation  problem.  Torques,  if  present,  contribute 
terms  of  0(R_i*)  in  the  drag.  Therefore,  Waklya's  (1965)  analyals  of  the 
drag  to  0(R  )  also  gives  the  sedimentation  velocity  to  at  least  that 
order. 

Waklya's  problem  is  equivalent  to  that  of  two  stationary  spheroids 

in  a  uniform  stream  with  the  stream  flowing  in  the  positive  x  direction 

(Figure  1).  We  start  by  deriving  the  method  of  reflections  solution  to 

Wakiya'a  resistance  problem,  but  without  any  simplifications  regarding 

relative  sizes,  or  spheroid  orientations  in  the  unifrom  stream  U*.  As 

in  Wakiya'a  work,  the  analysis  in  this  section  will  be  oarried  out  to 

two  reflections  so  that  the  drag  will  be  accurate  to  0(R~^).  The 

orientation  vector,  position  along  the  axis,  eccentricity  and  the 

distance  from  the  centroid  to  the  foci  of  each  spheroid  will  be  denoted 

by  d„*  *nd  c  ,  a-1,2. 
a  a  a 

In  an  earlier  note  (Kim  1984)  it  was  shown  that  the  Chwang-Wu 
(1974,1975)  representation  for  the  reflection  from  the  uniform  stream 
and  the  contribution  from  this  (zero-th)  reflection  to  the  drag  on 
spheroid  1  are: 


V*>  “  y"  "  U*-(mId1d1  ♦  0,(6-^)}. 
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The  Chwang-Wu  constants  a  and  Y  dtpsnd  only  on  ths  spheroid 
scesntricitiss  and  art  given  in  Table  1  and  I  is  the  Oseen  tensor.  The 
analogous  reflection  field,  y£  from  spheroid  2  can  be  obtained  from 
(2.1)  by  permuting  the  particle  Indices.  As  shown  in  the  refernces  on 
tfr~  method  of  reflections,  this  implies  that  the  contribution  from  the 
first  reflection  is 
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It  is  apparent  that  the  dependence  on  the  orientation  and  shape  of 
spheroids  1  and  2  comes  solely  from  the  tensor 


I(51-52)/(8w). 


(2.M) 


For  the  spherical  case,  this  tensor  is  known  to  polymer  kinetlcists  as 
the  Rotne-Prager-Yamakawa  tensor  (Rotne  and  Prager  1969,  Yamakawa  1970). 
It  is  interesting  to  note  that  Rotne  and  Prager  obtained  their  tensor 
using  a  variational  approach. 


The  leading  order  tern  in  the  contribution  fron  the  second 
reflection, 

cones  fron  the  nonopole  approxination, 

ij2  ~  36{a^^i 4  J(?r?2)/(8ir»,>  (2.5) 

The  dreg  on  spheroid  1  is  the  sun  of  these  contributions  and 
Waklya's  (1965)  solution  is  recovered  after  the  appropriate 
simplifications  in  the  geometry  and  notational  changes.  The 
contribution  fron  the  first  reflection  to  0(R  )  is  simply  a 


nonopole-nonopole  interaction. 

-20,(0^,^  ♦  et(«-d1d1)}*  I^-Xg)*^05  (2.6) 

To  oonvert  this  into  Waklya's  expression,  we  need 
F<0)  -  Ifivpc^ajd^  ♦  et(«-d2d2))«U*f  (2.6a) 

U"*(x2-x1 )  -  U*Rsint,  (2.6b) 

d^ »u  ■  d2«u  ■  o,  (2.6o) 


and  oro2-o,  e^eg-e. 

Equation  (2.6)  can  then  be  alnplifed  to 


and  the  x1,y1,z1  components  reduce  to  the  corresponding  terms  In 
Waklys's  (3.8). 

The  original  expressions  for  the  torque  can  also  be  recovered  via 
the  following  contributions  from  the  first  and  second  reflections. 
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In  the  next  section,  the  analogous  expressions  for  the 
sedimentation  (mobility)  problem  are  calculated  to  higher  order.  From 
here  on,  we  will  not  use  the  angles  t*  and  Instead,  the 
geometrical  dependence  will  be  represented  by  dot  products  between 
?1-x2,  d,  and  d2. 


3.  SEDIMENTATION  OF  TWO  SPHEROIDS 
3.1  General  Procedure 

The  procedure  for  calculating  the  sedimentation  velocities  is  a 

straightforward  generalization  of  that  employed  for  spheres.  The 

essential  modification  is  the  distribution  of  the  singularities  along 

the  axes  of  the  spheroids.  The  calculations  will  be  performed  up  to  and 

including  the  second  reflection  so  that  the  error  in  the  translational 

“7  —8 

and  rotational  velocities  will  be  0(R  )  and  0(R  )  respectively. 

The  translational  velocity  of  spheroid  1  is  obtained  by  summing  the 
contribution  from  the  reflections. 
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The  zero-th  reflection  contributes 

S}0*  “  ♦  al1  <^-d1d1 )} 

The  contribution  from  the  first  reflection  is 


with  the  incident  field 
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This  contribution  can  be  simplified  as  follows: 
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^th  g12  -  grg2  and  C  -  |g12|. 


The  second  reflection  contributes 


(3.6) 


with 


Yi2(x>-  (c|-q)n+(c|-5|)(^ffj.7«}vT(x-52)/(8*u)<152.  (3.7) 


When  we  Insert  (3.7)  into  (3.6)  end  use  the  expressions  for  the  Stokes 
dipole  end  octupole,  (3.6)  simplifies  to 
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The  stresslet  is  obtslned  vis  the  Psxen  lew  ss  (Kim  1984) 

S2ij  "  *  ?*,(d2id2j“  3*ij)(d2kd2r  ^kl5 

-  U«  <d21«jkd24  +  d2i6Jld2k  +  6ild2Jd2k  *  6ikd2Jd2l  ‘  i,d2id2Jd2kd2t) 

-  ^o„($ik6Jl  +  6il6Jk  +  6iJd2kd24  +  d2id2Jd2kd2l 

"  d2i6Jkd21  ’  d2i*Jtd2k  *  4ild2Jd2k  “  6ikd2Jd2l)  1X 

L°c2  (c|-6|)  {1  ♦  (cI-Cp^V)  e1kl(52)852 


2»yY  (d2lcJk4<l2l  ♦  d2JCikld215  *-ol  (c2"52)7x^1  (52^^52 * 


(3.9) 


The  expression  for  uj1^  is  exact  since  the  exact  y2  is  used  while 
is  accurate  only  to  0(R  8)  since  only  those  leading  terms  were  used  in 
(3.7).  However,  the  leading  error  term,  of  0(R*7),  comes  from  the  third 
reflection  whioh  was  neglected. 
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The  rotational  velocity  of  spheroid  1  also  follows  as  a  sum  of  the 


contribution  from  each  reflection. 
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The  first  reflection  contributes 
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Substitutions  for  y2,  its  rate-of- strain  field  g2  and  the  expression  for 
the  rotlet  eventually  lead  to 
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The  second  reflection  contributes 
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The  leading  order  term  in  y12  is  irrotational.  Therefore,  the  0(R  ■*) 
term  in  the  rotational  velocity  comes  solely  from  the  second  term  in 
(3.13).  It  should  be  noted  that  this  term  is  absent  for  spheres,  so  for 
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spheres,  u}  decays  more  rapidly,  l.e.,  as  0(R  r)  (see  Jeffrey  and 

Oniahi  0984)).  Finally,  (3*13)  can  be  reduced  to 


(3.14) 


Equations  (3.2),  (3.5),  (3.8),  (3.9),  (3.12)  and  (3*14)  allow  us  to 
oalculate  the  translational  and  rotational  velooitles  for  any  given 
orientations  of  the  spheroids.  Special  orientations  are  examined  In  the 
following  subsections. 

3.2  Sedimentation  of  two  vertically  oriented  spheroids  along  their  line 
of  centers 

The  simplest  geometry  consists  of  two  vertically  oriented  spheroids 
falling  along  their  line  of  centers  as  shown  In  Figure  2.  The 
orientation  and  separation  between  the  spheroids  remain  invariant  as  can 
be  deduced  from  symmetry  and  the  linearity  of  the  Stokes  equation.  Thus 
the  analysis  leads  to  a  straightforward  extension  of  Stlmson  and 
Jeffery’s  (1926)  results  for  spheres.  Therefore  this  problem  merely 
serves  as  a  benchmark  test  for  computational  methods. 


FIGURE  2.  Sedimentation  of  two  identical  spheroids  falling  along  their 
line  of  centers. 

The  sedimentation  velocity,  nondlmensionalized  by  the  terminal 
velocity  for  an  isolated,  vertically  oriented  spheroid,  is  plotted  in 
Figure  3  for  three  aspect  ratios.  The  solid  line  is  indistinguishable 
from  the  spherical  (Stimson- Jeffery)  solution  for  R/a  >2.1.  As  the 
aspect  ratio  increases,  the  hydrodynamic  interactions  become  weaker. 
The  contribution  from  reflections  beyond  those  calculated  by  Waklya 
become  significant  for  R/a  <  3*  The  results  also  agree  with  those 
obtained  by  Gluckman  et.  al.  (1971)  and  Liao  and  Krueger  (1980)  as 
shown  in  Table  2  of  the  appendix.  Their  1  factor,  the  drag 
nondlmensionalized  by  the  Stokes  drag  of  the  sphere  with  the  same 
cross-sectional  area,  has  been  successfully  reproduced. 
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SEDIMENTATION  VELOCITY 


Figure  4.  Sedimentation  of  two  horizontally  and  parallely  oriented 
spheroids . 

3.3  Sedimentation  of  two  horizontally  oriented  spheroids 
The  analysis  of  two  spheres  falling  perpendloular  to  their  line  of 
oenters  can  be  extended  to  the  spheroidal  ease.  However,  to' keep  the 
geometry  invariant,  the  spheroid  axes  must  be  perpendloular  to  both 
gravity  and  the  line  of  oenters. 

Figures  5  and  6  show  the  sedimentation  and  angular  velocities  for 

4 

the  same  aspect  ratios  as  before.  The  terminal  veloelty  of  an  isolated, 

horizontally  oriented  spheroid  was  used  to  scale  both  -11  and  -u  a.  In 

*  y 

both  figures  the  solid  line  is  indistinguishable  from  the  result  for 
spheres  (Goldman,  Cox  and  Brenner  1966 ) J  The  dependence  on  the  aspect 
ratio  is  qualitatively  similar  to  that  found  in  the  previous  subsection 
Finally,  it  should  be  noted  that  for  the  geometry  considered  in  this 
subseotlon,  the  0(R~^)  term  vanishes  in  equation  (3 .13)  for  the 
rotational  velocity. 


SEDIMENTATION  VELOCITY 


FIGURE  6.  Angular  velocity  vs.  center-to-center  separation  for  two 

spheroids  as  in  Figure  2.  Aspect  ratio  of  1.01,  - ; 

aspect  ratio  of  2,  - ;  aspect  ratio  of  10,  - — -. 
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3.4  Sedimentation  of  two  inclined  spheroid* 

The  results  in  the  preceding  subsections  were  qualitatively  similar  to 
that  found  for  spheres,  and  were  mainly  of  Interest  as  benchmarks  for 
the  computational  technique.  Zn  this  section,  we  turn  our  attention  to 
a  situation  where  the  results  differ  qualitatively  because  of  the 
evolution  of  the  particle  geometry.  Figure  7  shows  two  inclined 
spheroids  settling  with  their  axes  lying  in  a  common  vertical  plane. 
Hlrror  symmetry  has  been  introduced  to  reduce  the  number  of  parameters, 
but  the  algorithm  from  Subsection  3.1  can  handle  more  general 
situations.  At  all  times,  the  geometry  is  specified  by  the 
dimensionless  center- to-center  separation,  R/a  and  e,  the  azimuthal 
angle  for  dj. 
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Figure  7.  Mirror  symmetry  geometry  of  two  inclined  spheroids  with  their 
axes  in  a  vertical  plane 


The  successive  improvement  obtained  with  each  new  reflection  is 
shown  in  Table  1  of  the  appendix.  At  all  but  small  separations,  our  two 
reflection  solution  provides  accurate  answers.  For  spheres  the  exact 
result  is  available  (Goldman  et.  al.  1966)  and  we  see  that  even  at 
fairly  small  separations,  the  relative  error  is  under  lOf  because  of  the 
small  contributions  from  the  neglected  terms. 

The  evolution  of  the  geometry  is  caused  by  the  anisotropy  in  the 
mobility  tensors  and  the  rotation  of  the  spheroids  about  their 
respective  minor  axes.  Since  the  mobility  is  greater  in  the  axial  than 
in  the  transverse  direction,  an  inclined  spheroid  drifts  horizontally  as 
it  settles.  At  the  same  time,  the  spheroid  rotation  changes  the 
orientation  of  the  axis.  These  two  effects,  under  the  quasi-steady 
assumption,  are  governed  by  the  dimensionless  equations  (with  R/a 
rewritten  now  as  R) 

«  -  uy(R,6)  (3.15) 

and 

R  -  -2Ux(R,e)  (3.16) 

Figures  8  and  9  show  the  evolution  of  R  and  6  as  determined  by 
integrating  (3.15)  and  (3*16)  with  a  fourth  order  Runge-Kutta  routine. 
The  plots  Include  the  curve 

R  -  2(1  -  e*cos*e)1/2, 

for  contact  between  the  two  spheroids.  Figures  10  and  11  show  the 
corresponding  trajectories  of  the  centroid  of  spheroid  2  in  the  x-z 
plane. 


ORIENTATION  TRAJECTORIES 


ORIENTATION  TRAJECTORIES 


on  and  separation  for  two  spheroids 
.ling  as  in  Figure  7. 


FIGURE  9.  Evolution  of  orientati 
aspect  ration  =  2,  fal 
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PARTICLE  TRAJECTORIES 


If  the  orientation  trajectories  are  followed  fr on  e-0  (vertically 
oriented  epheroide)  and  all  allowed  valuea  for  R,  the  curves  in  Figures 
8  and  9  fall  into  two  groups,  depending  on  the  initial  value  of  R.  If  R 
exceeds  a  critical  value  at  6-0,  the  particles  eventually  and 
aonotonleally  separate  and  the  orientations  approach  asymptotically  a 
limiting  value  for  6  which  is  less  than  v/2  because  at  large 
separations,  g  goes  to  zero.  However,  for  initial  values  of  R  less  than 
the  critical  value,  the  rotational  motion  is  sufficiently  large  to  cause 

the  particles  to  rotate  beyond  the  horizontal  orientation.  Thereafter, 

« 

the  particles  drift  towards  eaoh  other  along  trajectories  which  are 
mirror  images  of  the  outward  trajectories.  The  separatrlx  which  starts 
at  the  critical  value  of  R  has  the  asymptote  e"-ir/2  (horizontal 
orientation) . 

The  asymptotic  behavior  at  large  R  can  be  obtained  by  using  only 
the  leading  term  in  «  and  that  in  U^.  Equations  (3*15)  and  (3.16)  can 
then  be  Integrated  analytically.  It  is  found  that  the  trajectories 
approach  the  limiting  orientation  6*  as 

R  1  -  (1  -  a1/a1)(cos26  -  cos2e")/(^ec1).  (3.17) 


The  influence  of  the  aspect  ratio  is  seen  by  comparing  Figure  8 

with  9  for  aspect  ratios  of  10  and  2  respectively.  For  slender 

particles,  the  periodic  trajectories  must  squeeze  through  a  very  narrow 

corridor  at  e-0.  As  the  aspect  ratio  is  reduced,  this  corridor  widens 
and  the  periodic  trajectories  become  more  like  the  straight  vertical 

lines  of  the  spherical  case. 

The  preceding  analysis  has  shown  that  R(t)  and  6(t)  for  the 
particle  geometry  of  this  subsection  are  either  periodic  or  represent 
single  encounters.  Accurate  calculation  of  the  separatrlx  required  the 
higher  reflections,  particularly  at  large  aspect  ratios. 
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APPENDIX 

The  numerical  convergence  with  each  additional  reflection  la  shown  in 
the  following  tables,  pie  geometry  is  as  in  Figure  7  and  the  velocities 
have  been  scaled  with  U  ,  the  sedimentation  velocity  of  an  isolated* 
vertically  oriented  spheroid. 


Table  1.  Inclined  Spheroids:  Effect  of  Suooesslve  Reflections 


Aspect  Ratio  *  1  and  R/a  -  3.0 

Uix/u‘ 

"iz^ 

.  “ly^ 

zero-th  reflection 
with  first  reflection 
with  second  reflection 
exact  solution 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

1.26852 

1.26852 

1.26680 

0.00000 

0.08333 

0.08333 

0.08178 

Aspect  Ratio  ■  2,  8  ■ 

0  and  R/a  • 

•  1.5 

U1x/U“ 

°1z/D‘ 

«1ya/U 

zero-th  reflection 
with  first  reflection 
with  second  reflection 

0.00000 

0.00000 

0.00000 

1.00000 

1.32692 

1.32886 

0.00000 

0.12419 

0.12457 

Aspect  Ratio  -  2,  6  • 

0.3*  and  R/a  »  2 

°1X/U" 

U1z/U" 

»iyW 

zero-th  reflection 
with  first  reflection 
with  second  reflection 

0.06033 

0.06033 

0.05296 

0.91697 

1.18429 

1.18259 

0.00000 

0.13897 

0.14073 

Aspect  Ratio  -  10,  9 

-  0  and  R/a 

-  0.7 

U1x/U" 

U1z/U* 

u»1y&/U 

zero-th  reflection 
with  first  reflection 
with  second  reflection 

0.00000 

0.00000 

0.00000 

1.00000 

1.28195 

1.28785 

0.00000 

0.07562 

0.06962 

Aspect  Ratio  -  10,  9 

-  0.3*  and  R/a  -  2 

"ix^ 

«iz/U- 

WiyS/U' 

zero-th  reflection 
with  first  reflection 
with  second  reflection 

0.14531 

0.14531 

0.14238 

0.80000 

0.91759 

0.91754 

0.00000 
6. 06406 
0.06416 

v>. 

$2 

mmmm 

'vV.sV-. 

r.v./V.V. 
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»  *  .  *\*^.  *  .  %*  '*  - 

i'’  ; 

t&k, 

i  >  v*  V  •/'  < 

Table  2.  Comparison  with  the  boundary  collocation  solution  of  Gluckman 
et.  al.  0971)  for  axisymaetrio  uniform  streaming.  X  is  the 
spheroidal  drag  divided  by  the  drag  on  a  sphere  with  the  same 
cross-sectional  area. 


Aspect  ratio  -  2 


R/a  Method  of  reflections  Collocation 

2  0.6485  0.8442 
4  0.9811  0.9812 
6  1.0458  1.0458 


Aspect  ratio  -  5 


R/a  Method  of  refleotions  Collocation 

2  1.3873  1.3700 
4  1.5675  1.5675 
6  1.6364  1.6364 
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NOTATION 


a  a* jor  seal-axis  of  spheroid, 

b  ainor  seal -axis  of  spheroid, 

c  distance  froa  center  to  foci, 

d  unit  vector  denoting  orientation  of  spheroid  axis 
e  eccentricity  of  the  spheroid, 

rate-of-straln  tensor, 
force  exerted  on  the  particle  by  the  fluid, 
gravitational  vector. 

Oseen  tensor. 

center  to  center  separation  between  two  spheroids 
stresslet  or  symmetric  part  of  the  stress-dipole, 
torque  exerted  on  the  particle  by  the  fluid, 
particle  translational  velocity., 
velocity, 
position  vector. 

constants  In  the  Chwang-Wu  singularity  solutions, 
constants  In  the  Chwang-Nu  singularity  solutions. 
Identity  tensor, 
alternating  tensor, 
angle  defined  In  Figure  7. 
viscosity. 

vector  denoting  position  on  the  spheroid  axis, 
stress  tensor. 

♦1#f2  angles  defined  In  Figure  1. 

♦  angle  defined  In  Figure  1. 

«  particle  angular  velocity. 


0  vortlclty. 

Q  vortlclty  tensor. 

Subscripts 

1,2  refers  to  spheroids  at  g|t 

1»J#M*«  Indices  used  In  the  Einstein  sumatlon  convention 
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(n)  denotes  association  with  the  n-th  reflection. 
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